Pattern Formation of Glioma Cells: Effects of Adhesion

We investigate clustering of malignant glioma cells. \emph{In vitro} experiments in collagen gels identified a cell line that formed clusters in a region of low cell density, whereas a very similar cell line (which lacks an important mutation) did not cluster significantly. We hypothesize that the mutation affects the strength of cell-cell adhesion. We investigate this effect in a new experiment, which follows the clustering dynamics of glioma cells on a surface. We interpret our results in terms of a stochastic model and identify two mechanisms of clustering. First, there is a critical value of the strength of adhesion; above the threshold, large clusters grow from a homogeneous suspension of cells; below it, the system remains homogeneous, similarly to the ordinary phase separation. Second, when cells form a cluster, we have evidence that they increase their proliferation rate. We have successfully reproduced the experimental findings and found that both mechanisms are crucial for cluster formation and growth.Comment: 6 pages, 6 figure

Predicting the Location of Glioma Recurrence After a Resection Surgery

International audienceWe propose a method for estimating the location of glioma recurrence after surgical resection. This method consists of a pipeline including the registration of images at different time points, the estimation of the tumor infiltration map, and the prediction of tumor regrowth using a reaction-diffusion model. A data set acquired on a patient with a low-grade glioma and post surgery MRIs is considered to evaluate the accuracy of the estimated recurrence locations found using our method. We observed good agreement in tumor volume prediction and qualitative matching in regrowth locations. Therefore, the proposed method seems adequate for modeling low-grade glioma recurrence. This tool could help clinicians anticipate tumor regrowth and better characterize the radiologically non-visible infiltrative extent of the tumor. Such information could pave the way for model-based personalization of treatment planning in a near future

Multiscale modelling of vascular tumour growth in 3D: the roles of domain size & boundary condition

We investigate a three-dimensional multiscale model of vascular tumour growth, which couples blood flow, angiogenesis, vascular remodelling, nutrient/growth factor transport, movement of, and interactions between, normal and tumour cells, and nutrient-dependent cell cycle dynamics within each cell. In particular, we determine how the domain size, aspect ratio and initial vascular network influence the tumour's growth dynamics and its long-time composition. We establish whether it is possible to extrapolate simulation results obtained for small domains to larger ones, by constructing a large simulation domain from a number of identical subdomains, each subsystem initially comprising two parallel parent vessels, with associated cells and diffusible substances. We find that the subsystem is not representative of the full domain and conclude that, for this initial vessel geometry, interactions between adjacent subsystems contribute to the overall growth dynamics. We then show that extrapolation of results from a small subdomain to a larger domain can only be made if the subdomain is sufficiently large and is initialised with a sufficiently complex vascular network. Motivated by these results, we perform simulations to investigate the tumour's response to therapy and show that the probability of tumour elimination in a larger domain can be extrapolated from simulation results on a smaller domain. Finally, we demonstrate how our model may be combined with experimental data, to predict the spatio-temporal evolution of a vascular tumour

The concentration of three anti-seizure medications in hair: the effects of hair color, controlling for dose and age

BACKGROUND: This paper assess the relationship between the quantity of three anti-seizure medications in hair and the color of the analyzed hair, while controlling for the effects of dose, dose duration, and patient age for 140 clinical patients undergoing anti-seizure therapy. Three drugs are assessed: carbamazepine (40 patients), valproic acid (40 patients), and phenytoin (60 patients). The relationship between hair assay results, hair color, dose, dose duration, and age is modeled using an analysis of covariance. The covariance model posits the hair assay results as the dependent variable, the hair color as the qualitative categorical independent variable, and dose, dose duration, and age as covariates. The null hypothesis assessed is that there is a no relationship between hair color and the quantity of analyte determined by hair assay such that darker colored hair will demonstrate higher concentrations of analyte than lighter colored hair. RESULTS: The analysis reveals that there is a significant relationship between dose and concentration for all hair color categories independent of the other covariates or the categorical independent variable. CONCLUSION: There does not appear to be any relationship between carbamazepine concentration and hair color. There is a weak relationship between hair color and valproic acid concentration, which the data suggest may be mediated by age. There is a significant, moderate relationship between phenytoin concentration and hair color such that darker colored hair has greater concentration values than lighter colored hair

Source localization of reaction-diffusion models for brain tumors

We propose a mathematically well-founded approach for locating the source (initial state) of density functions evolved within a nonlinear reaction-diffusion model. The reconstruction of the initial source is an ill-posed inverse problem since the solution is highly unstable with respect to measurement noise. To address this instability problem, we introduce a regularization procedure based on the nonlinear Landweber method for the stable determination of the source location. This amounts to solving a sequence of well-posed forward reaction-diffusion problems. The developed framework is general, and as a special instance we consider the problem of source localization of brain tumors. We show numerically that the source of the initial densities of tumor cells are reconstructed well on both imaging data consisting of simple and complex geometric structures

Determining and interpreting correlations in lipidomic networks found in glioblastoma cells

Background: Intelligent and multitiered quantitative analysis of biological systems rapidly evolves to a key technique in studying biomolecular cancer aspects. Newly emerging advances in both measurement as well as bio-inspired computational techniques have facilitated the development of lipidomics technologies and offer an excellent opportunity to understand regulation at the molecular level in many diseases. Results: We present computational approaches to study the response of glioblastoma U87 cells to gene- and chemo-therapy. To identify distinct biomarkers and differences in therapeutic outcomes, we develop a novel technique based on graph-clustering. This technique facilitates the exploration and visualization of co-regulations in glioblastoma lipid profiling data. We investigate the changes in the correlation networks for different therapies and study the success of novel gene therapies targeting aggressive glioblastoma. Conclusions: The novel computational paradigm provides unique “fingerprints” by revealing the intricate interactions at the lipidome level in glioblastoma U87 cells with induced apoptosis (programmed cell death) and thus opens a new window to biomedical frontiers. Background Glioblastoma are highly invasive brain tumors. Th

Tumor Angiogenesis and Vascular Patterning: A Mathematical Model

Understanding tumor induced angiogenesis is a challenging problem with important consequences for diagnosis and treatment of cancer. Recently, strong evidences suggest the dual role of endothelial cells on the migrating tips and on the proliferating body of blood vessels, in consonance with further events behind lumen formation and vascular patterning. In this paper we present a multi-scale phase-field model that combines the benefits of continuum physics description and the capability of tracking individual cells. The model allows us to discuss the role of the endothelial cells' chemotactic response and proliferation rate as key factors that tailor the neovascular network. Importantly, we also test the predictions of our theoretical model against relevant experimental approaches in mice that displayed distinctive vascular patterns. The model reproduces the in vivo patterns of newly formed vascular networks, providing quantitative and qualitative results for branch density and vessel diameter on the order of the ones measured experimentally in mouse retinas. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of different parameters in this process, hence underlining the necessary collaboration between mathematical modeling, in vivo imaging and molecular biology techniques to improve current diagnostic and therapeutic tools

A general reaction-diffusion model of acidity in cancer invasion

We model the metabolism and behaviour of a developing cancer tumour in the context of its microenvironment, with the aim of elucidating the consequences of altered energy metabolism. Of particular interest is the Warburg Effect, a widespread preference in tumours for cytosolic glycolysis rather than oxidative phosphorylation for glucose breakdown, as yet incompletely understood. We examine a candidate explanation for the prevalence of the Warburg Effect in tumours, the acid-mediated invasion hypothesis, by generalising a canonical non-linear reaction–diffusion model of acid-mediated tumour invasion to consider additional biological features of potential importance. We apply both numerical methods and a non-standard asymptotic analysis in a travelling wave framework to obtain an explicit understanding of the range of tumour behaviours produced by the model and how fundamental parameters govern the speed and shape of invading tumour waves. Comparison with conclusions drawn under the original system—a special case of our generalised system—allows us to comment on the structural stability and predictive power of the modelling framework

A mathematical modelling tool for predicting survival of individual patients following resection of glioblastoma: a proof of principle

The prediction of the outcome of individual patients with glioblastoma would be of great significance for monitoring responses to therapy. We hypothesise that, although a large number of genetic-metabolic abnormalities occur upstream, there are two ‘final common pathways' dominating glioblastoma growth – net rates of proliferation (ρ) and dispersal (D). These rates can be estimated from features of pretreatment MR images and can be applied in a mathematical model to predict tumour growth, impact of extent of tumour resection and patient survival. Only the pre-operative gadolinium-enhanced T1-weighted (T1-Gd) and T2-weighted (T2) volume data from 70 patients with previously untreated glioblastoma were used to derive a ratio D/ρ for each patient. We developed a ‘virtual control' for each patient with the same size tumour at the time of diagnosis, the same ratio of net invasion to proliferation (D/ρ) and the same extent of resection. The median durations of survival and the shapes of the survival curves of actual and ‘virtual' patients subjected to biopsy or subtotal resection (STR) superimpose exactly. For those actually receiving gross total resection (GTR), as shown by post-operative CT, the actual survival curve lies between the ‘virtual' results predicted for 100 and 125% resection of the T1-Gd volume. The concordance between predicted (virtual) and actual survivals suggests that the mathematical model is realistic enough to allow precise definition of the effectiveness of individualised treatments and their site(s) of action on proliferation (ρ) and/or dispersal (D) of the tumour cells without knowledge of any other clinical or pathological information